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PREFACE 


This  report  summarizes  the  work  completed  during  the  first  year  of  a  project  focused  on  im¬ 
proving  the  WBMOD  model  of  ionospheric  scintillation.  This  work  is  part  of  a  larger  efTort 
which  has  the  overall  objective  of  providing  the  USAF  Air  Weather  Service  with  an  im¬ 
proved  capability  for  providing  support  to  their  customers  who  are  impacted  by  ionospheric 
scintillation  effects. 

We  would  like  to  thank  Dr.  Santimay  Basu  of  the  USAF1  Phillips  Laboratory  (PL/GP1A) 
for  providing  the  scintillation  data  sets  from  Ascension  Island,  Huancayo,  and  Manila,  and 
for  his  guidance  and  suggestions  for  improving  the  model. 


IX 


1.  Introduction 


Many  modern  military  systems  used  for  communications,  command  and  control,  naviga¬ 
tion,  and  surveillance  depend  on  reliable  and  relatively  noise  free  transmission  of  radiowave 
signals  through  the  earth’s  ionosphere.  Small-scale  irregularities  in  the  ionospheric  density 
can  cause  severe  distortion,  known  as  radiowave  scintillation,  of  both  the  amplitude  and 
phase  of  these  signals.  A  basic  tool  used  in  estimating  these  effects  on  systems  is  a  computer 
program,  WBMOD,  based  on  a  single-scatter  phase-screen  propagation  model  and  a  number 
of  empirical  models  of  the  global  morphology  of  ionospheric  density  irregularities.  A  recent 
validation  of  the  WBMOD  model  showed  that  it  was  deficient  in  a  number  of  areas.  Much 
of  this  deficiency  could  be  traced  to  the  limited  coverage  of  the  data  used  in  developing 
the  model  (a  single  station  at  high  latitudes  and  two  stations  at  equatorial  latitudes).  The 
objective  of  this  study  is  to  investigate  and  implement  improvements  to  the  WBMOD  model 
based  on  more  extensive  data  sets  in  both  the  equatorial  and  high-  latitude  regimes. 

2.  Scintillation  Database  Development 

The  data  to  be  used  in  the  modeling  effort  come  from  two  sources.  The  first  is  data 
collected  from  the  Wideband  satellite  in  the  1970s  at  Ancon,  Peru,  and  Kwajalein  Atoll. 
The  second  is  data  collected  from  the  Marisat  satellite  from  Ascension  Island  and  Huancayo, 
Peru,  and  from  the  Fleetsat  satellite  from  Manila,  Philippines.  The  second  data  set  was 
provided  by  PL/GPIA.  Both  data  sets  were  used  in  a  validation  of  the  WBMOD  scintillation 
model,  reported  in  Secan  [1991].  Table  1  lists  the  coverage  of  these  data  sets.  The  column 
headed  “IPP”  is  the  geomagnetic  (apex)  latitude  of  the  350-km  ionospheric  penetration 
point  (IPP).  For  Ancon  and  Kwajalein,  this  is  defined  for  the  overhead  case. 


Station 

Satellite 

IPP 

Dates 

Points 

Ancon 

Wideband 

0.9 

1976-1979 

28,550 

Huancayo 

Marisat 

1.6 

1976-1984 

332,800 

Kwajalein 

Wideband 

4.1 

1976-1979 

30,420 

Manila 

Fleetsat 

5.8 

1981-1985 

Ascension  Island 

Marisat 

-16.5 

1980-1989 

Table  1:  Data  used  in  the  equatorial  model  development. 

2.1  Converting  SI(dB)  to  S4.  The  first  step  in  preparing  the  SI(dB)  data  provided  by  PL 
for  use  in  modeling  was  to  convert  it  to  the  equivalent  S4  value.  This  was  to  be  done  using 
an  algorithm  described  in  Secan  [1991]  (taken  from  Whitney  [1974])  based  on  an  empirically 
derived  relationship  between  the  57  (d/7)  index  and  the  Nakagami  m-parameter  (Figure  6  in 

1 


Whitney  [1974]).  It  is  not  clear  over  what  range  of  m  values  this  relationship  is  valid,  in 
particular  if  it  is  valid  for  low  values  of  m  (which  translates  into  high  values  of  S’*).  1  his 
algorithm  uses  the  following  equation  to  calculate  S 4  from  SI(dB): 

S<  =  0.0735S/(d/J)°757  (1  ) 

Our  concern  over  the  validity  of  this  algorithm  was  first  raised  in  an  analysis  of  data 
from  Ascension  Island  which  was  used  in  the  WHMOD  validation  study  ( Secan  [1991]).  I  he 
result  of  this  analysis  was  that  it  appeared  that  the  S4  values  derived  from  the  Ascension 
Island  SI(dI3)  values  were  roughly  20%  too  high  (see  the  Appendix  in  Secan  [1991]).  A 
solution  to  this  problem  was  approached  from  two  directions.  First,  empirical  factors  were 
derived  which,  when  multiplied  times  the  S4  values  generated  from  the  SI(dB)  data,  led  to 
distributions  of  S4  at  the  PL  stations  which  were  in  agreement  with  the  distributions  seen  in 
the  Wideband  data  sets.  Second,  a  set  of  digital  data  from  Ascension  Island  was  analyzed 
by  PL/GPIA  to  calculate  both  S4  and  SI(dB)  to  be  used  to  test  the  empirical  relationship 
generated  from  the  information  in  Whitney  [1974]. 

The  three  histogram  plots  in  Figure  1  illustrate  the  problem.  These  plots  show  the' 
distribution  of  S4  from  the  Kwajalein  and  Ancon  Wideband  data  sets  (top  and  bottoms  plots, 
respectively)  and  the  distribution  of  S4  calculated  from  SI(dB)  from  the  Ascension  Island 
MARISAT  L-band  data  set  (middle  plot).  The  data  used  in  generating  these  distributions 
spanned  the  time  period  from  one  hour  prior  to  local  sunset  to  eight  hours  after  sunset  . 
Note  that  the  distributions  for  the  two  Wideband  sets  peak  just  below  an  S4  value  of  1.0  (at 
0.93  for  both  sets)  caused  by  saturation  in  the  level  of  intensity  scintillation  due  to  Fresnel 
filtering.  The  similar  peak  in  the  Ascension  Island  data  is  at  a  value  of  1.10.  While  it 
is  possible  for  S4  to  exceed  1.0  in  cases  of  strong  multiple  scatter,  it  is  unlikely  that  the 
peak  in  the  distribution,  which  is  caused  by  the  saturation  effect,  would  occur  above  1.0. 
Table  2  summarizes  the  locations  of  the  peaks  in  the  various  S4  data  sets.  Note  that  the  two 
Huancayo  data  sets  peak  at  values  above  the  peaks  in  the  Wideband  data  sets,  and  that  the 
Manila  data  set  peaks  below  the  Wideband  value. 

It  is  interesting  to  note  that  if  we  regard  the  S4  value  of  0.93  from  both  directly-computed 
Wideband  VHF  data  sets  as  the  expected  S4  value  at  this  peak,  a  “corrective”  multiplier 
of  0.95  (±8%)  for  Huancayo  VHF  data  and  1.09  (±6%)  for  Manila  VHF  data  is  needed. 
Although  the  data  were  insufficient  for  a  true  64  peak  in  the  Wideband  L-band  data,  if 
we  used  the  Wideband  VHF'  peak  of  0.9d  also,  a  corrective  multiplier  of  0.97  is  needed  for 
the  Huancayo  L-band  data.  All  of  these  adjustments  are  quite  small,  and  are  essentially 
indistinguishable  from  no  correction,  given  the  gross  assumptions  used  in  developing  the 
value  of  the  corrective  multiplier. 

This  multiplier  for  the  Ascension  Island  data  set,  however,  is  0.84.  This  is  significantly 
different  from  1.0,  and  is  in  very  good  agreement  with  the  value  of  0.8  derived  from  a 
comparison  of  the  Ascension  Island  SI(dB)  values  with  published  samples  of  the  digital 
Ascension  Island  data  [Secan,  1991]. 
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Kwajalein  -  WIDEBAND  (VHF) 


Ancon  -  WIDEBAND  (VHF) 


0.0  0.5  1.0  1.5  2.0 

S4 


Data  limited  to  -lh  <  E-region  sunset  <  +8h 


Figure  i:  5. i  distributions  from  the  Kwajalein  and  Ancon  Wideband  VHF  data  sets  (top 
and  bottom  plots)  ami  derived  from  the  Ascension  Island  Marisat  L-band  SI(rlB)  data  set. 
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Station 

Frequency 

Satellite 

.S’ 4  method 

Peak  location 

Kwajalein 

VHF 

V  ideband 

direct 

0.93  (  +  /-  0.025) 

Ancon 

VHF 

Wideband 

direct 

0.93  (  +  /-  0.025) 

Huancayo 

VHF 

Marisat 

convert  from  SI 

0.98  (  +  /-  0.05) 

Manila 

VHF 

Fleetsa* 

convert  from  SI 

0.85  (+/-  0.025) 

Kwajalein 

L-  band 

Wideband 

direct 

insufficient  data 

Ancon 

L-band 

Wideband 

direct 

insufficient  data 

Huancayo 

L-band 

Marisat 

convert  from  SI 

0.95  (+/  0.025) 

Ascension 

L-band 

Marisat 

convert  from  SI 

1.10  (+/■  0.025) 

Table  2:  Location  of  peaks  in  S4  distribution  plots. 

In  the  second  approach,  personnel  at  PL/GPIA  analyzed  a  small  sample  of  digital  inten¬ 
sity  data  collected  at  Ascension  Island  and  calculated  both  S4  and  SI{dD)  from  the  same 
15 -minute  samples.  These  data  were  fit  with  an  equation  in  the  same  form  as  Equation  [1] 
relating  S4  to  SI(dB),  which  resulted  in  the  following  relationship: 

=  0.0764S/(dfl)°72  (2) 

Figure  2  shows  the  relationship  between  SI(dB)  and  S4  from  the  original  Whitney  algorithm 
(dotted  line),  from  the  fit  to  the  Ascension  Island  digital  data  (solid  line),  and  from  the 
Whitney  algorithm  modified  with  the  multiplier  of  0.84  derived  from  the  S4  distributions 
(dashed  line).  The  agreement  between  the  results  obtained  from  the  digital  data  and  that 
from  the  comparison  of  the  54  distributions  is  excellent. 

As  a  result  of  these  analyses,  we  have  decided  to  use  the  Whitney  algorithm  modified 
by  the  multipliers  derived  from  the  locations  of  the  peaks  in  the  S4  distributions  to  convert 
SI(dD)  to  Si .  We  have  not  attempted  to  determine  why  the  unmodified  Whitney  algorithm 
did  not  work  correctly,  nor  why  the  adjustment  factors  are  different  at  the  different  stations 
(and  for  different  frequencies).  It  is  possible  that  these  factors  represent  calibration  factors 
which  are  required  to  bring  the  various  measurements  into  relative  agreement  with  one 
another.  It  is  enough  for  our  purposes  that  the  adjustment  agrees  with  an  algorithm  derived 
directly  from  the  data  in  the  one  case  for  which  the  data  are  available  (Ascension  L-band), 
and  that  we  apply  it  consistently  across  our  data  sets. 

2.2  Database  Design.  The  equatorial  database  files  designed  to  facilitate  the  WBMOI) 
model  validation  conducted  in  1991  are  large  and  cumbersome,  and  they  do  not  lend  them 
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Comparison  of  SI-to-S4  Algorithms 


SI(dB) 


Figure  2:  Comparison  of  algorithms  for  converting  SI[dB)  to  54- 

selves  to  the  types  of  studies  envisioned  for  the  modeling  process.  We  have  redefined  the 
layout  of  the  equatorial  database  tailored  to  the  types  of  modeling  we  expect  to  do  and  the 
parameters  we  expect  to  require  for  modeling. 

The  validation  database  stored  each  observation  as  an  independent  entity.  This  leads 
to  a  somewhat  disjoint  database  and  quite  a  bit  of  data  duplication  (both  implicitly  and 
explicitly).  We  reformatted  the  databases  provided  by  Phillips  Lab  (Huancayo.  Manila, 
Ascension  Island)  so  that  the  data  are  stored  in  data  blocks  which  will  include  24  hours  of 
data  to  include  the  SI(dB)  observations  as  well  as  parameters  to  be  used  in  the  modeling. 
Figure  3  is  a  listing  of  the  design  of  a  data  structure  for  this  new  database.  A  single  data 
structure  contains  th"  SI(dB)  values  for  a  24-hour  period  starting  with  the  first  value  after 
local  solar  noon  on  a  given  day.  Definition  of  each  element  is  as  follows: 


INTSECO  Date/time  of  first  SJ(dB)  value  (seconds) 

SI (96)  SI(dB)  values  for  the  24-hour  period 

ELTO  Time  since  E-layer  sunset  of  first  SI(dB)  value  (hours) 

TPSSO  Time  since  E/F2-laycr  sunset  of  first  Sl{dB)  value  (hours) 

DOY  Day  of  the  yea*.  (GMT)  of  first  SI(dB)  value 

PHITRM  Angle  between  sunset  terminator  and  the  geomagnetic  meridian 

KP(9)  Kv  values  coveiing  the  24-hour  period 

KP.SS  A'r  value  at  F/F2  sunset 
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STRUCTURE  /DATA  RECORD/ 


INTEGER 

INTSECO 

REAL 

SI (96) 

REAL 

ELTO 

REAL 

TPSSO 

REAL 

DOT 

REAL 

PHITRM 

REAL 

KP  (9) 

REAL 

KP_SS 

REAL 

SUM6_KP_SS 

REAL 

DST (25) 

REAL 

dst_min 

REAL 

slt_dst_min 

REAL 

AVE_DST_DAT 

REAL 

AV£_DST_PSS 

REAL 

AVE_DST_PM 

REAL 

SSN 

REAL 

SSN_5 

REAL 

SSN_90 

REAL 

FIO 

REAL 

F10_5 

REAL 

F10_90 

REAL 

SPARES (2) 

END  STRUCTURE 


igure  3:  Database  structure  for  data  from  Ascension  Island,  Huancayo,  and  Manila. 
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SUM6_KP_SS 

DST(25) 

DST.MIN 

SLT_DST_MIN 

AVE_DST_DAY 

AVE_DST_PSS 

AVE_DST_PM 

SSN 

SSN.5 

SSN_90 

F10 

F10.5 

F 10.90 

SPARES (2) 


52  Kv  for  six-hour  period  prior  to  E/F2- layer  sunset 
I),t  values  covering  the  24- hour  period 
Minimum  Dat  observed  during  the  24-hour  period 
Local  solar  time  (LS  I  )  of  Dat  minimum  (hours) 
Average  Dat  between  local  noon  and  E/F2  sunset 
Average  Dal  between  E/F2  sunset  and  local  midnight 
Average  Dat  between  local  midnight  and  0600  LST 
Sunspot  number  (SSN) 

Five-day  mean  SSN 

90-day  mean  SSN 

10.7  cm  solar  radio  flux  (F10) 

Five-day  mean  F10 

90- day  mean  F10 

Unused  space  for  future  use 


Note  that  the  Wideband  data  sets  from  Ancon  and  Kwajalein  were  kept  in  the  form  used 
in  the  validation.  These  data  are  from  individual  satellite  passes,  and  the  database  format 
used  in  the  validation  is  better  for  this  type  of  data  than  the  new  format.  We  will  also  retain 
the  validation  format  for  the  high-latitude  modeling  for  the  same  reason. 

The  various  geomagnetic-activity  parameters  included  in  this  list,  such  as  the  various 
average  Dat  values,  are  based  on  our  literature  review  to  date.  They  were  selected  to  provide 
parameters  which  would  reflect  observed  behavior,  such  ac  in  Aarons  (1991),  and  theoretical 
modeling,  such  as  in  Kelley  and  Maruyama  [1992],  reported  in  the  literature. 

3.  Irregularity  Models  Development 

3.1  Model  for  Iog(C*L)  PDF.  From  its  inception,  the  WBMOD  model  has  dealt  strictly 
with  the  mean  of  the  CkL  distribution.  While  this  treatment  is  reasonable  for  mid  and 
high  latitudes  where  the  log (Ci,L)  distribution  is  very  nearly  normal  (i.e.,  the  distribution  is 
well  represented  by  a  Gaussian  curve),  recent  studies  have  shown  that  it  is  not  adequate  for 
representing  conditions  in  the  equatorial  post-sunset  regime  (Secan  and  Fremouur[1988j).  As 
shown  in  Figure  4,  the  distribution  of  log(C*L)  in  this  regime  is  rather  more  bimodal  than 
normal,  which  leads  to  the  situation  where  the  mean  of  the  distribution  is  not  representative 
of  the  distribution.  In  addition,  characterizing  only  the  mean  of  the  distribution  does  not 
allow  for  providing  model  users  with  statistics  which  may  be  of  use  to  them  such  as  generation 
of  954,*-percentile  values  for  worst-case  estimates  or  the  percent  of  the  time  that  scintillation 
levels  will  exceed  a  user-defined  threshold. 

In  the  previously  cited  study,  it  was  found  that  the  log(CVZ,)  distributions  observed  at  ail 
latitude  regimes  (equatorial,  mid  latitude,  auroral,  and  polar  cap)  could  be  modeled  using  a 
linear  combination  of  Gaussian  functions.  In  all  but  the  equatorial-night  case,  only  a  single 
function  was  required.  In  the  equatorial-night  case  (as  in  Figure  4),  the  distribution  could 
be  adequately  modeled  by  a  three-component  PDF. 
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Figure  4:  Distribution  of  log((’t/,)  during  the  post-sunset  period  at  Kwajalein  and  Ancon 
from  the  Wideband  experiment  (VHF). 

The  functional  form  of  the  model  to  be  used  for  the  log (CtL)  probability  distribution 
function  (PDF)  will  be  a  composite  function  consisting  of  twogaussian  peaks  joined  together 
by  a  “top  hat'’  function  which  has  a  fixed  value  between  the  peaks  and  is  zero  outside  this 
range.  The  functional  form  is  given  I)}’ 

PDF  =  ]£(?/,  -  ?/,;)exp  -  ^  +r/„  (3) 

>h  if  Hi  <  log (CkL)  <  H7 

1 b  (4) 

0  elsewhere 

where  (//,,  cr;, i  =  1,2)  define  the  !og(C'jt/>)  values,  the  half-widths,  and  the  distribution- 
peak  values  for  the  two  peaks,  and  is  the  value  of  the  “top  hat”  segment  of  the  model. 
After  the  various  regional  models  are  used  to  determine  values  for  these  seven  parameters,  the 
resulting  PDF  is  normalized  so  that  it  integrates  to  1.0  ( i.e CDF[ log(6’jt/>)  — >  co]  — ►  1.0). 
Figure  5  is  an  example  of  a  PDF  and  the  corresponding  CDF  constructed  using  this  new 
model  formulation.  It  can  be  compared  with  the  actual  distributions  shown  in  Figure  4. 

The  modeling  now  lo  be  completed  will  describe  the  variation  of  these  seven  parameters 
as  functions  of  location,  time,  and  geophysical  conditions. 
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Figure  5:  Probability  distribution  function  (PDF)  (top  plot)  and  cumulative  distribution 
function  (CDF)  (bottom  plot)  for  an  equatorial  station.  These  were  generated  using  the 
new  three-component  PDF  model  (see  text)  with  the  parameters  indicated  in  the  side  box. 
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3.2  Modeling  Philosophy.  As  we  began  to  analyze  the  data  to  extract  the  variations  in 
the  seven  parameters  described  above,  it  became  evident  that  we  would  need  to  alter  the 
way  in  which  we  have  developed  the  parameters  for  the  WBMOD  C*  L  model  in  the  past 
(see  Robins  et  al.  [1986]  and  Secan  (1989}  for  disc  ussions  of  previous  procedures).  Several 
factors  have  motivated  this  change: 

1.  We  are  now  modeling  the  full  distribution  of  log (CkL)  rather  than  simply  the  mean  of 
the  distribution.  Thus,  there  are  far  more  interrelated  parameters  to  deal  with. 

2.  Intensity  scintillation,  which  makes  up  the  majority  of  the  data  to  be  used,  saturates 
at  a  level  which  is  frequency-dependent.  In  addition,  the  hand-scaled  SI(dB)  data 
have  artificial  “quantization”  effects  near  the  noise  floor  (the  location  of  which  is  also 
frequency  dependent).  These  two  effects  make  it  impossible  to  extract  the  parameters 
required  directly  from  £4  (or  log5<)  PDFs. 

3.  Due  to  the  range-limited  nature  of  the  intensity  scintillation  data  (i.e.,  noise-floor  at 
the  low  end  and  saturation  at  the  high),  it  can  be  difficult  to  determine  whether  an 
increase  in  the  number  of  observations  near  the  top  of  the  £4  distribution  is  caused 
by  an  increase  in  the  population  in  the  plume  distribution  (an  increase  in  »/2)  or  an 
increase  in  the  log(Cjt^)  value  at  the  peak  (an  increase  in  p2). 


After  much  trial  and  error,  we  have  come  up  with  the  following  procedures  and  “rules 
of  thumb”  for  the  modeling  effort: 

1.  The  half-widths  of  the  two  gaussians,  <J\  and  02,  will  be  set  at  values  derived  from 
examination  of  the  log (C'kL)  PDFs  generated  from  the  WIDEBAND  04,  data  sets. 
Since  the  phase  scintillation  does  not  saturate  at  high  levels  of  C*b,  these  data  provide 
the  best  picture  of  the  shape  of  the  higher  post-sunset  gaussian  peak  associated  with 
plume  structures.  In  addition,  since  the  WIDEBAND  data  were  machine  processed 
rather  than  hand  scaled  and  are  available  at  a  lower  frequency  (137  MHz),  they  also 
provide  a  better  picture  of  the  shape  of  the  lower  gaussian  peak. 

2.  The  log(CfcL)  of  the  lower  gaussian  peak,  fx  1,  will  also  be  taken  from  the  WIDEBAND 
VHF  data  set  for  the  reasons  stated  in  the  last  sentence  of  the  previous  paragraph. 

3.  In  order  to  determine  whether  to  increase  or  /12  (as  in  item  #3  above),  we  will 
follow  the  simple  rule  of  thumb  that  if  we  believe  that  the  variation  being  studied 
should  result  in  an  increase  in  AN,  we  will  model  the  variation  by  changing  p2.  If 
we  believe  that  the  variation  is  related  to  the  probability  of  the  existence  of  a  plume 
structure,  we  will  model  the  variation  by  changing  rj2.  We  have  tried  to  hold  to  this 
rule  as  much  as  possible,  although  it  is  sometimes  difficult  to  decide  which  of  these  two 
cases  are  being  dealt  with. 
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4.  Since  the  r;,  values  of  the  resulting  PDF  are  to  be  normalized  so  that  the  PDF  integrates 
to  unity,  we  will  set  the  value  of  i) t  to  a  fixed  value  and  model  changes  in  the  relative 
sizes  of  the  upper  and  lower  populations  by  altering  the  value  of  r/2-  In  order  to  simplify 
the  modeling  effort,  we  have  also  decided  to  set  the  value  of  the  “top  hat”  segment  of 
the  model,  »/3 ,  to  a  fixed  percentage  of  q2. 

5.  The  modeling  will  be  conducted  by  comparing  contour  plots  of  the  percent  of  the  time 
that  exceeds  two  threshold  values  (0.5  and  0.8  for  all  but  the  Huancayo  L-band 
data;  0.15  and  0.25  for  Huancayo  L-band)  as  generated  from  the  various  data  sets  and 
as  generated  from  the  model.  In  as  many  cases  as  possible,  the  contour  plots  will  be 
a  function  of  GMT  along  the  x-axis  and  day  of  the  year  along  the  y-axis.  Figure  6 
is  an  example  of  such  a  plot  from  the  Manila  data  set.  The  constants  of  the  models 
developed  to  describe  the  variations  of  the  various  PDF  parameters  will  be  iteratively 
adjusted  to  provide  visual  best  fits  to  the  entire  data  set. 


As  stated  in  item  #3,  we  have  tried  to  follow  these  rules  and  procedures  as  much  as 
possible,  although  some  deviations  have  been  required.  Usually,  when  we  have  found  that  a 
particular  rule  or  procedure  was  not  working  well,  it  was  either  modified  or  replaced.  The 
modeling  procedure  has  definitely  grown  and  changed  in  the  course  of  the  modeling  effort. 

3.3  Latitudinal  Variation.  The  algorithm  for  the  latitudinal  variation  of  CyL  across  the 
equatorial  region  in  the  current  WBMOD  model  is  a  very  simple  one,  which  does  not  allow 
for  higher  CyL  values  near  the  crest  of  the  anomaly  region  than  at  the  equator.  In  order  to 
allow  for  this  sort  of  variation  with  latitude,  wc  will  need  to  replace  this  algorithm  with  one 
which  will  permit  this  sort  of  behavior.  The  current  model  specifies  CyL  in  the  mid-latitude 
region  and  at  the  geomagnetic  equator,  then  uses  a  simple  error  function  to  transition  from 
one  to  the  other  across  a  specified  boundary  latitude  with  a  specified  transition  width.  The 
new  algorithm  will  be  similar,  but  will  add  gaussian- shaped  peaks  at  the  poleward  edges  of 
the  equatorial  region  (north  and  south)  which  are  defined  by  a  new  set  of  model  constants. 


The  form  of  the  new  latitudinal  algorithm  is  as  follows: 


=  Go  +  (Cc  —  Co)  exp 


A„  1  —  Af 


H 

7c  =  { 


0  |  Att  |>  A c 

Ceo  j  Aa  j<  Ac 


Aa  >0 
Aa  <  0 
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f  WN  Aa  >  Ac 
Wc=  {  WB  |  Att  |<  Ac 
[Ws  Aa  <  -Ac 

where  A„  is  the  geomagnetic  (modified  apex)  latitude;  Ceo  is  the  equatorial  value  of  the 
parameter  varied;  Ccjv  and  Ccs  are  tire  maximum  values  of  the  parameter  varied  in  the 
northern  and  southern  anomaly  crest  regions;  Ac  is  the  geomagnetic  latitude  of  the  anomaly 
maxima  (always  positive);  and  W/v,  IVB,  and  Ws  are  the  e- folding  widths  of  the  gaussian 
features  north  of  +Ac,  equatorward  of  ±Ac,  and  south  of  —Ac-  All  of  these  parameters 
may  vary  as  functions  of  other  variables  such  as  time  of  day  or  season,  although  it  is  likely 
that  at  least  a  few  (such  as  the  transition  widths)  will  either  be  fixed  or  will  vary  as  simple 
functions  of  the  other  parameters.  Figure  7  shows  three  examples  of  the  latitude  variations 
possible  with  this  formulation. 

We  have  placed  the  following  initial  constraints  on  the  modeling  effort: 


1.  Only  the  log(CfcZ,)  value  of  the  upper  peak  in  the  PDF,  will  be  varied  with  latitude. 
This  is  tantamount  to  assuming  that  the  increase  in  scintillation  activity  levels  observed 
near  the  anomaly  crest  (at  Ascension  Island,  in  the  present  data  set)  is  due  to  increased 
levels  of  AN.  Since  this  has  been  the  primary  argument  for  the  increase  (i.e.,  there 
is  more  plasma  (N)  available  at  the  anomaly  crest  than  on  the  equator,  therefore  AN 
can  be  higher  at  the  crest),  we  felt  that  this  was  a  reasonable  first  assumption.  [As  we 
will  see  in  a  later  section,  this  had  to  be  modified  slightly.] 

2.  The  maximum  values  in  /i2  will  occur  at  apex  latitude  ±16°  with  a  latitude  half- width 
of  5°  on  each  side  of  the  peak.  The  location  of  the  peak  was  selected  to  coincide  with 
the  latitude  of  Ascension  Island,  the  only  data  in  the  set  near  the  crest,  and  the  width 
was  selected  to  reduce  to  the  equatorial  value  by  the  latitude  range  covered  by  the 
Wideband  data  sets,  which  showed  no  latitude  variation  in  CfcL. 

3.  The  latitude  variation  will  be  assumed  identical  in  all  longitude  sectors.  We  have  some 
reservations  about  this  assumption  based  on  studies  that  have  shown  that  the  higher 
scintillation  levels  observed  at  Ascension  may  have  a  longitude  component  as  well  as  a 
latitude  one  (see  Basu  el  al.  [1976]  and  Maruyama  and  Matuura  [1984]),  but  as  we  have 
no  other  scintillation  data  set  near  the  crest  from  another  longitude  sector,  we  must 
assume  (for  now)  that  the  differences  in  observed  scintillation  levels  at  Huancayo  and 
Ascension  Island  are  representative  of  the  latitude  variation  at  all  longitude  sectors. 
Note,  however,  that  this  is  definitely  an  area  that  needs  further  study  with  a  more 
extensive  data  set. 

4.  We  will  assume  (at  least  initially)  that  the  maximum  CkL  observed  at  Ascension 
corresponds  to  the  log(Cfc/j)  calculated  from  an  SI(dB)  value  of  30.0  (36.0),  which 
is  at  the  top  end  of  t>»  SI(dB)  distribution  from  that  station.  Due  to  the  saturation 
effect  found  in  intensity  scintillation,  this  could  well  be  low.  However,  we  have  no 
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I'igure  7:  Samples  of  the  now  equatorial  C\L  model  formulation.  The  model  parameters 
used  in  each  are  defined  in  the  adjacent  boxes.  See  the  text  for  a  description  of  the  model 
and  model  parameters. 
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evidence  to  set  log {C^L)  higher.  This  would  require  SI{dB)  observations  from  a  higher 
frequency,  or  phase  scintillation  observations  at  L-band  or  higher. 

Thus,  we  are  specifying  only  the  values  for  at  the  dip  equator  and  at  the  anomaly 
crest.  We  will  use  the  maximum  \og(CkL)  from  the  Ascension  Island  set  as  the  value  at 
the  crest,  and  a  value  of  33.5  at  the  dip  equator.  This  last  value  was  initially  taken  from 
the  average  value  of  the  log(CVL)  peak  in  the  two  Wideband  data  sets  [Secan,  1989].  and 
was  later  found  to  be  in  fair  agreement  with  a  value  extracted  from  the  Huancayo  UHF  and 
L-band  data  sets. 

3.4  Diurnal  Variation.  The  first  issue  to  resolve  in  modeling  the  diurnal  variation,  a 
variation  dominated  by  the  rapid  onset  of  intense  scintillation  just  after  local  sunset,  was 
to  select  the  time  metric  to  use  in  the  model.  The  current  WBMOD  equatorial  CkL  model 
uses  the  time  since  sunset  has  occurred  in  the  E-region  at  both  ends  of  the  field  line  passing 
through  the  ionospheric  penetration  point  (IPP).  This  time  metric,  which  we  will  denote 
ERSS  (for  E-region  sunset),  was  based  on  the  hypothesis  by  Tsunoda  and  Wittwer  (  Tsunoda 
[1985])  that  the  field-line  integrated  Pedersen  conductivity  was  a  major  controlling  factor 
in  the  onset  and  growth  of  the  instabilities  that  develop  into  post-sunset  plumes.  This 
parameter  is  largely  determined  by  the  electron  density  at  the  E-region  ends  of  a  particular 
field  line,  and  the  conductivity  drops  to  nearly  zero  shortly  after  both  E-region  ends  of  the 
field  line  have  passed  through  the  sunset  transition. 

The  development  of  the  diurnal  variation  in  the  WBMOD  equatorial  C*L  model  was 
severely  limited  by  the  fact  that  the  data  available  for  the  modeling  work  (Wideband  data 
from  Ancon  and  Kwajalein)  were  limited  to  a  two- hour  period  centered  at  4.0  ERSS  at.  Ancon 
and  a  ninety-minute  period  centered  at  4.5  ERSS  at  Kwajalein  (see  Figure  14  in  Fremouw 
and  Robins  [1985]).  The  data  sets  provided  by  PL/GPIA  cover  the  full  24  hours,  which  will 
allow  us  to  improve  this  section  of  the  model.  In  fact,  an  initial  improvement  was  made 
based  on  the  Ascension  Island  data  during  the  model  validation  (see  page  4  in  Secan  [1991]). 

The  plots  in  Figures  8  through  12  show  the  diurnal  behavior  of  the  SI  data  from  Ascension 
Island  (Marisat,  L-band  and  VHF),  Manila  (Fleetsat,  VHF),  and  Huancayo  (Marisat,  VHF 
and  L-band),  as  a  function  of  GMT  and  day  of  the  year.  These  figures  are  contour  plots 
of  the  percent  of  time  that  SI  exceeded  a  particular  threshold  (3  dB  for  Huancayo  L-band, 
10  dB  for  the  other  sets)  as  a  function  of  GMT  and  the  day  of  the  year  for  an  eight-hour 
time  segment  beginning  near  local  sunset.  The  heavy  solid  line  on  each  plot  is  the  time  of 
E-region  sunset  as  defined  above  (i.e.,  ERSS  =  0.0),  and  the  heavy  dashed  line  is  the  time  of 
sunset  in  the  F2  region.  (This  last  time  will  be  denoted  F2RSS,  so  the  heavy  dashed  line  is 
at  F2RSS  =  0.0.)  These  times  are  calculated  at  locations  to  be  defined  later.  It  is  apparent 
from  the  data  that  the  onset  of  scintillation  lags  the  F.RSS=G  time,  and  appears  to  track 
fairly  well  with  a  combination  of  the  two  sunset  times. 

The  decision  to  include  the  F2RSS=0  line  in  these  plots  was  based  on  recent  work  by 
Maruyama  [1988]  and  Kelley  and  Maruyama  [1992],  in  which  a  physical  model  for  the  growth 
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Figure  8:  Diurnal  behavior  of  intensity  scintillation  from  the  Ascension  Island  Marisat  (VHF) 
data  set.  Contours  of  the  percent  of  the  time  that  SI  exceeded  10  dB  are  plotted  as  a  function 
of  GMT  and  the  day  of  the  year  The  heavy  solid  line  shows  the  variation  of  ERSS  =  0  with 
the  day  of  the  year,  and  the  heavy  dashed  line  shows  the  variation  of  \  =  107°. 
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Figure  9:  Diurnal  behavior  of  intensity  scintillation  from  the  Ascension  Island  Marisat  (L- 
band)  data  set.  Contours  of  the  percent  of  the  time  that  SI  exceeded  10  dB  are  plotted  as  a 
function  of  GMT  and  the  day  of  the  year.  The.  heavy  solid  line  shows  the  variation  of  ERSS 
=  0  with  the  day  of  the  year,  and  the  heavy  dashed  line  shows  the  variation  of  \  —  107°. 
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Figure  10:  Diurnal  behavior  of  intensity  scintillation  from  the  Manila  Fleetsat  (VHF)  data 
set.  Contours  of  the  percent  of  the  time  that  SI  exceeded  10  dB  are  plotted  as  a  function  of 
GMT  and  the  day  of  the  year.  The  heavy  solid  line  shows  the  variation  of  ERSS  =  0  with 
the  day  of  the  year,  and  the  heavy  dashed  line  shows  the  variation  of  \  =  107°. 
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Figure  11:  Diurnal  behavior  of  intensity  scintillation  from  the  Huancayo  Marisat  (VHF)  data 
set.  Contours  of  the  percent  of  the  time  that  SI  exceeded  10  dB  are  plotted  as  a  function  of 
GMT  and  the  day  of  the  year.  The  heavy  solid  line  shows  the  variation  of  ERSS  =  0  with 
the  day  of  the  year,  and  the  heavy  dashed  line  shows  the  variation  of  \  =  107°. 
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Figure  12:  Diurnal  behavior  of  intensity  scintillation  from  the  Huancayo  Marisat  (VHF)  data 
set.  Contours  of  the  percent,  of  the  time  that  SI  exceeded  10  dB  are  plotted  as  a  function  of 
GMT  and  the  day  of  the  year.  The  heavy  solid  line  shows  the  variation  of  ER.SS  =  0  with 
the  day  of  the  year,  and  the  heavy  dashed  line  shows  the  variation  of  \  =  107°. 


of  post-sunset  plume  structures  was  developed.  This  model  indicates  that  photoionization 
in  the  F2  layer  is  an  important  “brake"  on  the  growth  process,  and  that  the  growth  rate 
increases  as  photoionization  in  the  F‘2  layer  decreases  at  sunset.  Thus,  as  our  initial  time 
metric,  we  have  chosen  a  sunset  time  for  onset  of  plume-structure  development  which  is 
a  combination  of  the  E-  and  F2- region  sunsets  described  earlier.  This  time,  denoted  tc  is 
defined  as  the  time  since  the  later  of  ERSS=0  or  F2RSS=0.  In  referring  to  Figures  8  through 
12,  te  is  F2RSS  on  those  days  when  the  dashed  line  is  to  the  right  of  the  solid  line  and  is 
ERSS  when  the  relative  positions  are  reversed.  For  Ascension  Island,  F2RSS  controls  the 
process  from  roughly  the  autumnal  equinox  to  the  vernal  equinox,  and  ERSS  controls  for 
the  rest  of  the  year.  At  the  other  two  stations,  the  F2RSS  controls  for  the  entire  year. 

In  our  initial  modeling,  we  used  our  standard  method  of  specifying  conditions  at  an 
ionospheric  penetration  point  (IPP)  for  calculating  the  various  sunset  times.  The  ERSS  and 
F2RSS  times  calculated  for  Huancayo  using  this  method  resulted  in  the  scintillation  onset 
preceeding  the  sunset  times  by  as  much  as  20  to  30  minutes.  (Note  that  the  time  lines 
plotted  on  Figures  II  and  12  were  calculated  using  the  new  procedure.)  This  early  onset 
was  a  concern  for  a  while,  but  it  was  eventually  resolved  in  principle  if  not  in  detail.  In  a 
recent  paper,  Mendillo  et  al.  [1992]  showed  two  sets  of  east- west  scans  from  the  Altair  radar 
located  at  Kwajalein  on  two  nights  in  August  1988.  These  scans  showed  several  examples  of 
the  plume  structures  that  are  the  source  of  the  most  severe  post-sunset  scintillation  events. 
In  particular,  the  scan  from  07:57  UT  on  14  August  (their  plate  la)  shows  an  intense  plume 
structure  located  at  the  eastern  boundary  of  the  scan,  and  the  scan  from  08:30  UT  on  15 
August  (plate  lb)  shows  a  plume  structure  located  just  cast  of  the  radar.  The  scans  just 
prior  to  both  of  these  (9  minutes  prior  on  14  August,  35  minutes  on  15  August)  were  devoid 
of  plume  signatures.  There  are  two  observations  that  can  be  made  about  these  plumes  that 
are  germane  to  the  question  at  hand  (timing  at  Huancayo): 

1.  Both  plumes  extend  well  above  the  sunset  terminator,  defined  as  a  line  passing  from 
the  center  of  the  sun  tangent  to  the  earth’s  surface.  A  consequence  of  this  is  that  a 
receiver  “looking”  eastward  at  a  low  elevation  angle  (such  as  the  Huancayo  Marisat 
geometry)  can  encounter  irregularities  in  the  topside  ionosphere  prior  to  the  terminator 
crossing  the  line  of  sight  at  ionospheric  altitudes. 

2.  In  the  non-plume  scans,  the  local  time  at  the  location  where  the  plume  develops  in 
the  next  scan  is  past  the  E-region  sunset  time  but  prior  to  the  F2  sunset  time.  In  the 
plume  scans,  the  local  time  at  the  plume  location  is  after  the  F2  sunset  time.  This 
is  in  accord  with  the  hypothesis  discussed  earlier  that  the  plumes  would  develop  after 
both  the  E  and  F2  sunset  times. 

This  provides  additional  information  that  may  allow  us  to  explain  and  model  the  variation 
in  onset  times  at  the  three  stations. 

The  F2  sunset  time  preceded  the  plume  scan  on  15  August  by  only  one  minute.  This  is 
not  enough  time  for  the  plume  to  have  developed  fully  to  a  height  of  over  800  km  (assuming 
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a  start  at  a  height  of  400  kin,  this  would  require  a  vertical  rise  velocity  of  almost  7  km/s). 
It  would  appear  that  the  time  that  photochemistry  in  the  F2  layer  no  longer  can  suppress 
the  growth  of  the  instability  requires  redefinition.  If  we  use  the  nine-minute  spacing  between 
the  two  scans  on  14  August  as  a  bound  and  parameterize  the  F2-layer  control  by  the  local 
solar  zenith  angle,  then  we  find  that  a  solar  zenith  angle  of  about  107°  corresponds  to  the 
earliest  time  that  the  plume  could  have  started  developing  in  the  14  August  case. 

Modeling  the  onset  time  will  require  consideration  of  the  vertical  extent  of  the  plumes 
rather  than  continue  to  work  solely  with  geometries  calculated  for  the  phase-screen  pene¬ 
tration  point.  The  geometry  at  the  three  stations  (Ascension,  Huancayo,  and  Manila)  are 
well  suited  to  investigating  this  issue,  as  the  propagation  paths  are  quite  different  at  each 
station.  At  Ascension,  the  azimuth  angle  of  the  raypath  is  nearly  along  the  geomagnetic 
meridian  at  a  high  elevation  angle;  at  Huancayo,  the  azimuth  is  eastward  (looking  towards 
the  approaching  sunset  terminator)  at  a  low  elevation  angle;  and  at  Manila,  the  azimuth  is 
westward  (looking  away  from  the  terminator),  also  at  a  relatively  low  elevation  angle. 

The  model  we  have  settled  on  for  the  zero  point  of  the  time  metric  to  be  used  in  the 
equatorial  region  will  be  the  time  that  the  E-region  points  along  the  field  line  passing  through 
the  equatorial  control  point  (definition  to  follow)  are  in  darkness  after  sunset  and  the  solar 
zenith  angle  at  the  latitude/longitude  of  the  control  point  is  >  107°.  This  time  metric  will 
be  denoted  te.  The  equatorial  control  point  is  determined  as  follows: 


1.  Calculate  the  UT  of  sunset  ( lt  =  0)  at  the  iatitude/iongitude  where  the  line  of  sight 
from  the  receiver  to  the  satellite  crosses  altitudes  defined  as  the  base  and  top  of  plume¬ 
generated  irregularities. 

2.  The  equatorial  control  point  is  the  one  of  these  two  locations  with  the  earliest  UT 
sunset  time. 

3.  The  field  line  along  which  the  E-layer  sunset  times  are  calculated  is  that  which  passes 
through  the  equatorial  control  point  at  the  base  altitude  of  the  plume  (as  in  (1)  above). 


Thus,  we  are  now  defining  the  time  in  terms  of  the  intersection  of  the  line  of  sight  with 
a  vertical  ‘‘wall”  of  irregularities  located  along  the  sunset  terminator  rather  than  with  the 
intersection  with  a  horizontal  phase  screen.  We  have  selected  altitudes  of  350  km  and  800  km 
for  the  base  and  top  of  the  plume  “wall”  based  on  qualitative  comparisons  with  the  Manila 
and  Huancayo  data  sets,  respectively. 

We  decided  to  retain  the  basic  functional  form  of  the  diurnal  variation  used  in  the  present 
version  of  WBMOD:  a  steep  rise  just  after  sunset  as  defined  above  using  a  gaussian  function, 
a  short-duration  flat  top  peak,  and  a  more  gradual  decrease  throughout  the  night  using  a 
second  gaussian  function.  We  will  dconte  the  two  gaussian  segments  the  sunset  and  sunrise 
segments.  The  actual  form  used  is  as  follows: 
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exe[-(^)2]  if 

/(*«)  =--  '  1.0  if  U<te<ttT  (6) 

.  exp  [-  i^Y]  if  *«  >  ^ 

where  <e  is  the  time  past  sunset,  and  VV„  are  the  peak  time  and  half-width  of  the  post¬ 
sunset  gaussian,  t„  is  the  time  at  which  the  “sunrise”  segment  of  the  variation  starts  its 
decrease,  and  WtT  is  the  half-width  in  the  sunrise  segment.  The  diurnal  variation  for  any 
particular  parameter  is  then  given  by 

X  ~  xday  +  ( xday  xnight)f{te)  (7) 

where  Xday  is  the  daytime  value  and  xnight  is  the  maximum  nighttime  value. 

We  found  in  modeling  that,  while  the  observed  variations  could  be  best  fit  by  using  the 
variation  modeled  in  Equation  (6]  for  the  Oi  and  r/,  parameters,  the  observed  rapid  onset 
of  severe  scintillation  that  is  occasionally  observed  at  all  stations  could  not  be  adequately 
modeled  unless  the  h?  parameter  (the  log(6’*L)  value  of  the  plume- associated  peak  in  the 
PDF)  was  set  to  its  maximum  value  immediately  upon  crossing  the  sunset  line.  In  order 
to  model  this,  we  have  set  the  diurnal  variation  function  for  the  two  Hi  parameters  to  the 
following: 


'  1.0 

if  te  < 

/,(*.)  =  '  10 

if  t'<te<  t,T 

(8) 

l  “P  [  te)  ] 

if  te  >  t„ 

This  function  is  then  used  in  place  of  f(tc)  in  Equation  (7]  for  the  diurnal  variation  of  the 
Hi  parameters. 

In  fitting  the  model  to  the  observations,  which  was  done  using  the  procedure  outlined 
in  Section  3.2,  it  was  found  that  while  the  tss,  W„ ,  and  t,r  parameters  could  be  set  as 
simple  constants,  the  W,r  parameter  needed  to  vary  as  a  function  of  te.  We  found  that  the 
scintillation  levels  dropped  off  more  rapidly  in  the  early  (pre-midnight)  evening  after  the 
peak  than  later  in  the  early  morning  hours.  We  found  that  we  were  unable  to  use  a  simple 
gaussian  with  a  constant  half-width  parameter.  In  order  to  keep  the  model  as  simple  as 
possible,  W,T  is  modeled  to  vary  linearly  with  <e  as  follows: 
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(9) 


W„  =  W.*  + 


t C 

dW.r 


In  addition,  we  found  that  this  variation  was  different  at  the  crest  (Ascension  Island)  than 
near  the  equator  (Huancayo  and  Manila)  so  there  are  different  W,r 0  and  dW,r  values  in  the 
equatorial  and  anomaly  crest  sections  of  the  model. 

The  set  of  parameters  for  the  diurnal  variation  is  as  follows: 


t.. 

1.50* 

t.r 

1.75* 

w„ 

0.80* 

o 

6.00* 

Equatorial 

2.00* 

Anomaly  Crest 

dW„ 

2.00* 

Equatorial 

1.50* 

Anomaly  Crest 

[Note:  The  overall  results  will  be  shown  in  Section  3.7.) 

3.5  Longitudinal  Variation.  It  was  found  that  in  order  to  reconcile  the  four  near¬ 
equator  data  sets  (Manila,  5.8°;  Kwajalein,  4.1°;  Ancon,  0.9°;  Huancayo,  1.6°)  a  longitu¬ 
dinal  variation  needed  to  be  included  in  the  k>g(C*L)  PDF.  This  has  been  included  in  the 
model  by  applying  a  longitudinal  modulation  on  the  rj2  parameter  such  that  the  value  of 
this  parameter  in  the  vicinity  of  Manila  is  a  factor  of  0.70  times  the  value  in  the  Kwa- 
jalein/Huancayo/Ancon/Ascension  sector.  The  very  non-uniform  longitudinal  coverage  of 
the  data  set  made  it  impossible  to  rigorously  define  this  variation,  so  it  has  been  imple¬ 
mented  in  the  following  simple  form: 


Vt, night  ~  Vmax  [0-70  +  0.30/(<£„)] 


(10) 


/(*»)  = 


«p[  (‘Sr1)! 

if  <t>„  <  <t> l 

1.0 

if  4>\  <  tf>v  <  <(> 2 

(ID 

«»l  (‘5?)1 

if  (j>„  >  <j>2 

where  <f>„  is  a  longitude  parameter  rotated  so  that  <f>n  =  0  occurs  at  a  geographic  longitude  of 
100°;  4>\  and  4>i  are  40°  and  320°,  respectively  (corresponding  to  140°  and  60°  E.  longitude); 
and  A <f>  is  set  to  10°.  These  values  were  set  in  order  to  move  the  transitions  away  from  the 
locations  where  we  had  data,  and  were  based  on  review  of  the  results  documented  in  Basu 
et  al.  [1976]  and  Maruyama  and  Matuura  (1984)  from  in-situ  observations. 
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3.6  Seasonal  Variation.  Since  we  are  now  modeling  the  full  log(G\ L)  distribution,  we  no 
longer  need  an  explicit  model  of  the  probability  of  the  occurrence  of  plumes  in  the  post¬ 
sunset  period.  Thus,  we  will  return  to  the  formalism  used  in  WBMOD  model  version  10111 
for  the  seasonal  variation  of  the  parameters  that  describe  the  log (CkL)  PDF.  Following  the 
modeling  philosophy  described  in  Section  ‘,1.2,  the  seasonal  behavior  is  modeled  via  alteration 
of  the  magnitude  of  the  upper  (plume)  peak  in  the  log((7*/,)  PDF'  (i.e.,  the  92  parameter). 
This  follows  from  the  assumption  that  the  seasonal  behavior  is  primarily  a  change  in  the 
probability  that  plumes  will  occur. 

As  in  version  10G1,  the  metric  used  in  modeling  the  seasonal  variation  is  the  angle 
between  the  geomagnetic  meridian  and  the  sunset  terminator,  4>term ,  with  the  maximum 
activity  expected  when  the  two  planes  defined  by  the  meridian  and  the  terminator  are  aligned 
=  0).  The  rationale  for  this  is  presented  in  Robins  et.  al  [1986).  The  functional  form 
of  the  seasonal  modulation  of  the  772  parameter  is  as  follows: 


f(<t>urm)  -  exp 


(12) 


where  W,  is  one  of  two  constants  --  one  for  local  summer  (15.0)  and  one  for  local  winter 
(12.0).  The  value  of  <t>uT7n  is  calculated  from  the  magnetic  declination  angle  and  the  angle 
to  the  terminator  (calculated  from  the  location  of  the  “true”  sun)  at  the  phase  screen  pen¬ 
etration  point.  The  determination  of  whether  to  use  the  summer  or  winter  value  of  W,  is 
also  made  at  that  location. 

3.7  Results  of  Diurnal/Seasonal/Longitudinal  Modeling.  The  results  of  the  modeling 
described  in  Sections  3.3  through  3.6  are  shown  in  Figures  13-22.  These  figures  are  in  pairs, 
with  the  first  figure  in  each  pair  showing  the  diurnal  and  seasonal  variation  of  the  percent 
of  time  that  S\  exceeds  the  thresholds  stated  in  Section  3.2  from  each  of  the  five  data 
sets  (Ascension  Island  L-band  and  VHF;  Iluancayo  VHF  and  L-band;  Manila  VHF),  and 
the  second  showing  the  same  percentages  as  calculated  from  the  model.  The  data  used  in 
generating  the  first  pair  of  plots  were  limited  to  times  of  low  geomagnetic  activity  ( Kp  <  2o), 
but  no  limits  were  placed  on  SSN. 

The  data  from  Ancon  and  Kwajalein  were  too  limited  in  local  time  to  make  comparable 
figures  for  those  two  data  sets.  In  comparing  Ancon  to  Huancayo,  we  found  (much  to  our 
relief)  that  when  the  Huancayo  VHF’  data  were  limited  to  the  times  and  conditions  covered 
by  the  Ancon  VHF  data,  the  variations  seen  in  the  two  data  sets  were  very  similar.  Thus, 
we  assume  that  if  the  model  behavior  is  correct  for  Huancayo,  it  is  probably  quite  good  for 
Ancon  as  well.  Kwajalein,  however,  is  another  story. 

Figure  23  shows  the  model  results  for  near-overhead  geometry  at  Kwajalein.  (The  two 
vertical  dotted  lines  indicate  the  GMT  range  covered  by  the  Kwajalein  data  set.)  This 
behavior  is  somewhat  similar  to  Manila,  but  with  more  scintillation  activity  in  the  local 
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Figure  14:  (a)  Model  results  for  Ascension  Island  (I, -band)  corresponding  to  the  data  shown 
in  13(a).  (b)  Model  results  for  Ascension  Island  (I, -hand)  corresponding  to  the  data  shown 
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Figure  16:  (a)  Model  results  for  Ascension  Island  (VHP)  corresponding  to  the  data  shown 
in  15(a).  (b)  Model  results  for  Ascension  Island  (VHP)  corresponding  to  the  data  shown  in 
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Figure  18:  (a)  Model  results  for  Iluancayo  (VHF)  corresponding  to  the  data  shown  in  !7(a) 
(|>)  Model  results  for  Iluancayo  (VHF)  corresponding  to  the  data  shown  in  17(b). 
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summer.  Figure  24  shows  a  comparison  of  the  post-sunset  data  from  Kwajalein  with  data 
from  Manila  limited  to  the  same  time-past-sunset  interval  as  the  Kwajalein  set  (roughly  0330 
to  0530  <e).  The  summer  “fill-in”  is  seen  to  be  far  greater  at  Kwajalein  than  at  Manila,  and 
greater  as  well  than  called  for  in  the  model.  This  was  also  noted  as  a  major  difference  between 
Kwajalein  and  Ancon  in  previous  modeling  [Secnn,  1989].  We  were  unable  to  reconcile  this 
feature  at  Kwajalein  with  the  observed  behaviors  at  the  other  stations  in  the  database, 
although  as  can  be  seen  in  Figure  23  the  model  does  tend  to  call  for  a  more  noticeable  “fill” 
in  the  summer  than  at  the  other  stations. 

In  obtaining  the  fits  as  shown  in  these  figures,  it  was  often  necessary  to  degrade  the 
quality  of  the  fit  at  one  station  (or  season,  or  time)  in  order  to  fit.  the  observed  behavior 
at  the  other  stations  (seasons,  times).  A  major  problem  of  this  type  can  be  noted  in  the 
peak  percentages  at  the  two  seasonal  maxima  at  the  various  stations.  In  each  case,  one  of 
the  maxima  is  higher  than  the  other.  Unfortunately,  there  is  no  systematic  rule  that  can  be 
extracted  from  the  data  in  hand.  The  higher  maximum  occurs  near  the  vernal  (post-summer) 
equinox  at  Huancayo  (and  Ancon),  the  autumnal  (pre-summer)  equinox  at  Ascension  Island, 
and  the  autumnal  (pre-winter)  equinox  at  Manila.  The  Kwajalein  data  show  a  broad  peak 
throughout  the  summer.  While  it  would  be  simple  to  add  a  seasonal  asymmetry  to  the 
occurrence  peaks  in  the  model  (requiring  only  a  minor  alteration  to  Equation  [12]),  the 
lack  of  a  systematic  behavior  has  forced  us  to  exclude  this  variation  and  simply  produce 
symmetric  peaks  with  peak  occurrence  values  midway  between  the  observed  values. 

3.8  Geomagnetic  Activity  Effects. 

WBMOD  has  never  included  effects  of  geomagnetic  activity  on  scintillation  levels  in  the 
equatorial  region.  The  reason  is  that  there  has  been  little  consensus  on  what  the  effects  are, 
aside  from  indications  that  increased  levels  of  activity  tend  to  surpress  post-sunset  scintilla¬ 
tion  and  may  enhance  post-midnight  scintillations.  Our  plan  of  attack  for  investigating  the 
possibility  of  adding  these  effects  is  as  follows: 

1.  Review  the  existing  literature  discussing  (a)  the  observational  evidence  and  (b)  current 
theories  of  the  effects  of  geomagnetic  activity  on  equatorial  scintillation. 

2.  Determine  which  parameters  to  use  in  partitioning  the  scintillation  data  and  in  the 
final  modeling.  For  example,  the  minimum  D,t  observed  during  the  previous  day,  Kp 
at  sunset,  £  Kp  over  the  previous  twelve  hours,  etc.  This  determination  will  be  based 
on  the  results  of  the  literature  review. 

3.  Conduct  both  statistical  and  case  studies  from  the  equatorial  database  of  the  behavior 
oi  scintillation  as  functions  of  the  various  geomagnetic  parameters. 

4.  If  a  systematic  behavior  is  found,  develop  an  empirical  model  for  the  effect  found. 

3.8.1  Literature  Review.  This  section  contains  the  results  of  a  preliminary  literature 
search  for  prior  research  work  on  the  subject  of  geomagnetic  control,  or  influence,  on  oqua- 
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tonal  scintillation.  It  must  be  noted  that  the  articles  cited  in  this  report  do  not  represent  a 
comprehensive  list  of  all  publications,  but  are  rather  a  selection  of  the  voluminous  papers  on 
the  subject.  Since  much  of  the  available  literature  examines  geomagnetic  effects  on  Equa¬ 
torial  Spread  F  (ESF),  and  the  ESF  phenomenon  is  intimately  related  to  equatorial  radio 
scintillation,  the  literature  search  will  consider  geomagnetic  control  or  influence  on  ESF'  as 
well  as  radio  scintillation.  There  is  considerable  disagreement  on  the  proper  measure  of  the 
level  of  geomagnetic  activity,  but  the  two  most  commonly  utilized  in  the  recent  literature 
seem  to  be  hp  and  Dlt ,  whereas  much  of  the  early  work  used  sunspot  activity  as  the  measure 
of  geomagnetic  activity.  As  background,  some  papers  demonstrating  an  understanding  of 
the  basic  processes  producing  scintillation  will  be  included,  even  where  direct  geomagnetic 
influence  on  scintillation  is  not  a  topic.  Early  work  on  equatorial  scintillation  often  used  radio 
star  sources  to  study  scintillation  phenomena,  [Wright  et  al.,  195G;  hoster,  1958;  hosier  and 
Wright,  I960].  These  works  focused  on  radio  scintillation  measurements  rather  than  on  ESF 
measurements  and  revealed  the  close  correlation  between  ESF  and  equatorial  scintillation, 
as  well  as  the  first  indications  of  a  negative  correlation  between  equatorial  scintillation  and 
increased  geomagnetic  activity. 

Satellite  scintillation  work  improved  the  studies,  and  early  work  strengthened  the  sugges¬ 
tion  that  high  geomagnetic  activity  often  has  an  inhibitory  efTect  on  equatorial  scintillation. 
Koster  [1972]  used  the  Kp  index  to  show  that  for  lower  values  of  Kp  (53  hp  <  30),  there  is 
no  clear  correlation  with  scintillation;  however,  for  large  values  of  hp  {^Zhp  >  30)  there  is 
a  strong  negative  correlation  with  scintillation.  Mullen  [1973]  compared  ATS-3  scintillation 
data  (SI)  and  Dsl  data,  and  showed  evidence  that  scintillation  occurrence  in  the  June  solstice 
increased  during  periods  of  magnetic  activity,  but  that  the  general  efTect  of  magnetic  activity 
appears  to  be  a  suppression  of  seasonal  scintillation  occurrence,  with  an  annual  distribution 
tending  toward  uniformity. 

Present-day  work  on  modeling  of  equatorial  scintillations  requires  a  clear  understanding 
of  the  basic  physical  mechanisms  at  work.  Perhaps  the  onset  of  this  understanding  was  when 
strong  evidence  was  presented  by  Woodman  and  La//oz[1976]  showing  that  ESF1  conditions 
and  radio  scintillation  effects  were  due  to  plume-shaped  electron  density  irregularities  or 
bubbles.  This  study  was  based  on  data  from  the  Jicamarca  radar  data  analysis.  This 
understanding  of  the  basic  mechanism  of  radio  scintillation  began  to  lead  to  much  better 
modeling  of  the  details  of  the  mechanisms  involved. 

Basu  and  helley  [1979]  reviewed  the  equatorial  scintillation  observations  and  their  rela¬ 
tionships  to  the  theoretical  development  of  plumes  and  bubbles,  but  did  not  include  magnetic 
activity  effects  upon  scintillation.  In  addition,  Ossakow  [1981]  reviewed  much  of  the  the¬ 
ory  and  the  numerical  simulation  results  for  plume  structure  development,  which  strongly 
increased  understanding  of  the  basic  mechanisms  of  plume  formation  due  to  the  Rayleigh- 
Taylor  instability  mechanism,  though  he  also  did  not  address  the  geomagnetic  activity  in¬ 
fluences. 

Some  early  observational  work  was  complicated  by  the  differences  between  the  “range” 
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and  “frequency”  types  of  equatorial  spread  F.  It  was  shown  in  Rastogi  (1980)  that  the  re¬ 
lationship  between  range  spread  F  and  sunspot  activity  is  positive  for  pre-midnight  hours, 
and  strongly  negative  for  post-rnidnight  sectors.  The  frequency  spread  F,  however,  showed 
relatively  little  sunspot  activity  effects. 

Aarons  et  al.  (1980)  looked  at  both  the  seasonal  and  geomagnetic  control  of  equato¬ 
rial  scintillations  in  two  longitudinal  sectors.  Although  the  authors  did  not  delve  into  the 
underlying  mechanisms,  they  found  that,  in  general,  pre-midnight  scintillation  is  inhibited 
by  magnetic  activity  (based  on  Iip)  and  that  during  the  post-midnight,  period  increased 
magnetic  activity  increases  scintillation  activity  under  moderate  solar  flux  conditions.  A 
review  of  the  status  of  ionospheric  scintillations  was  given  by  Aarons  (1982).  Although  the 
empirically  derived  magnetic  activity  effects  upon  scintillation  were  described,  the  actual 
mechanisms  of  control  were  still  undefined.  Since  the  prime  driver  of  the  Rayleigh-Taylor 
mechanism  is  the  ExB  drift  after  sunset  forcing  the  layer  to  rise  and  become  unstable,  it  is 
clear  that  modeling  the  geomagnetic  effects  on  scintillations  will  require  an  understanding 
of  the  equatorial  electric  fields  involved.  Fejer  (1986)  presented  a  summary  of  the  equatorial 
electric  fields  associated  with  magnetic  disturbances.  He  found  that  the  zonal  equatorial 
electric  fields  often  show  large  and  sudden  perturbations  during  geomagneticaily  disturbed 
periods,  whereas  the  F-region  vertical  electric  fields  do  not  change  significantly  with  magnetic 
activity.  He  also  detailed  that  to  first  order,  perturbation  electric  fields  are  westward  during 
the  day  and  eastward  at  night,  and  that  southward  Bt  (component  of  the  Interplanetary 
Magnetic  Field)  changes  seem  less  effective  in  disturbing  the  equatorial  electric  fields. 

Earle  and  Kelley  (1987)  have  shown  that  for  Iip  >  3,  fluctuations  in  the  magnetospheric 
electric  fields  with  periods  less  than  about  six  hours  penetrate  virtually  unattenuated  to 
the  equatorial  ionosphere,  while  those  with  longer  periods  are  shielded  out.  A  clear  and 
concise  account  of  the  high  latitude  effects  on  the  equatorial  electrodynamics  is  given  in 
Kelley  (1989). 

Aarons  (1991)  introduced  a  view  that  the  ring  current  (as  determined  by  the  peak  excur¬ 
sion  of  D,t)  during  magnetospheric  storms  plays  a  leading  role  by  either  directly  or  indirectly 
establishing  the  conditions  necessary  for  equatorial  F  layer  irregularity  generation  or  inhi¬ 
bition.  The  view  that  the  relationship  between  the  local  time  and  the  time  of  maximum 
ring-current  energy  is  of  primary  importance  to  the  development  of  plumes  is  emphasized. 
The  theory  says  that  if  the  maximum  ring-current  energy  (as  shown  by  l),t)  occurs  during 
the  midnight  to  postmidnight  time  period,  irregularities  are  generated.  If  the  maximum  l),t 
occurs  in  the  early  afternoon,  then  irregularities  are  inhibited.  If,  however,  the  maximum 
occurs  around  sunset  or  shortly  after  sunset,  then  there  is  no  effect  on  the  generation  of  ir¬ 
regularities  that  night.  The  behavior  of  scintillation  during  a  number  of  individual  magnetic 
storm  events  were  presented  to  support  this  view.  The  weaknesses  of  this  paper  are  that 
only  individual  events  were  presented,  no  statistical  analysis  of  a  large  data  set  being  given, 
and  that  the  connection  between  the  maximum  ring-current  energy  and  the  F-layer  electric 
fields  is  neither  dearly  presented  nor  justified  by  observational  data. 
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Recently  a  paper  by  Kelley  and  Maruyama  [1992]  presented  a  diagnostic  model  for  equa¬ 
torial  spread  F  detailing  the  effect  of  magnetic  activity.  This  work  is  based  upon  the  key 
observation  that  for  short  periods  (less  than  six  hours),  high-latitude  electric  fields  penetrate 
to  low  latitudes  (see  previous  section).  The  basic  model  used  was  developed  in  a  previous 
paper  by  Maruyama  [1988].  In  this  paper,  by  means  of  a  numerical  model,  they  examined 
the  electric  field  effects  on  the  generation  of  ESF  in  the  post-midnight  hours. 

The  above  model  was  designed  to  run  using  average  zonal  electric  field  patterns  published 
by  Fejer  et  al.  [1979],  based  on  season  and  solar-cycle  variations.  A  gaussian-like  perturba¬ 
tion  was  added  to  simulate  an  eastward  electric  field  penetration  from  high  latitudes  to  the 
average  zonal  electric  field.  This  perturbation  was  added  in  different  model  runs  at  different 
times  during  the  night,  and  it  was  found  that  an  increase  in  the  instability  growth  rate  (7) 
was  not  just  due  to  the  destabilizing  influence  of  an  eastward  component  of  the  electric  field. 
The  altitude  change  of  the  layer  and  the  corresponding  decrease  in  the  ion-neutral  collision 
frequency  (1/)  (which  increased  the  value  of  the  gju  term)  was  an  important  factor.  The 
g/v  term  dependends  on  the  ratio  of  the  gravitational  field  to  the  collision  frequency.  The 
model  runs  showed  that,  in  general,  the  later  in  the  evening  the  disturbance  occurs,  the 
more  effective  it  is  in  enhancing  instability  growth.  As  would  be  expected,  the  perturbation 
amplitude  also  plays  an  important  part  in  determining  the  growth  rate.  Thus,  the  model 
shows  that  a  reversal  of  the  electric  field  from  zonally  westward  to  zonally  eastward,  while 
the  main  contributor  to  and  initiator  of  the  generation  mechanism,  does  not  have  absolute 
and  direct  control.  Indeed,  the  reversal  is  sometimes  not  sufficient  by  itself  to  trigger  intense 
ESF.  The  altitude  change  of  the  layer  and  the  corresponding  decrease  in  ion-neutral  collision 
frequency  has  a  large  efTect  in  the  model,  and  often  is  the  crucial  parameter  when  determined 
by  the  magnitude  and  duration  of  the  event. 

The  Kelley  and  Maruyama  [1992]  study  proceeds  to  examine  individual  cases,  which 
are  quantitatively  very  well  modeled.  Electric  field  data  from  Jicamarca  was  used  so  that 
comparison  to  actual  plume  development  could  be  performed.  The  caveat  was  presented 
that  the  observational  data  was  limited  since  the  flux  tubes  involved  were  so  long  that 
no  single  station  could  model  all  necessary  parameters.  The  transequatorial  thermospheric 
wind  effects  discussed  in  the  Maruyama  [1988]  paper  were  also  not  included  in  the  Kelley 
and  Maruyama  [1992]  paper  since  no  observational  data  existed  for  the  days  in  question. 
Although  no  statistical  database  was  used  for  analysis,  the  model  does  explain  the  previous 
observations:  that  spread  F  which  develops  post-midnight  does  occur  more  often  during 
periods  of  magnetic  disturbance. 

The  status  of  modeling  the  efTect  of  geomagnetic  activity  upon  equatorial  spread  F  and 
equatorial  radio  scintillations  is  improving.  It  is  clear,  however,  that  there  is  a  tremendous 
amount  of  work  to  be  done  yet  upon  this  subject,  and  the  time  when  such  models  can  be 
used  as  actual  predictors  of  scintillation  is  still  far  away. 

3.8.2  Model  Development.  In  developing  the  way  in  which  WBMOD  will  reflect 
the  effects  of  geomagnetic  activity  on  equatorial  scintillation,  we  held  to  the  philosophy 
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that  WBMOD  is  a  climatological  (i.e.,  statistical)  rather  than  a  phenomenological  model. 
In  other  words,  the  model  will  attempt  to  reproduce  the  behavior  seen  in  the  occurrence 
statistics  rather  than  attempting  to  model  it  in  terms  of  a  set  of  “canonical”  storm-time 
behaviors.  We  have  also  not  attempted  to  include  many  of  the  details  of  the  behaviors 
we  have  found  in  the  data,  in  most  part  because  it  was  almost  impossible  to  find  a  set  of 
behaviors  common  to  all  stations  beyond  that  which  we  have  decided  to  implement. 

In  keeping  with  the  modeling  philosophy  outlined  in  Section  3.2,  we  will  model  the  effects 
of  geomagnetic  activity  on  equatorial  scintillation  in  terms  of  its  effect  on  the  probability 
that  plumes  will  form  and  grow.  This  will  be  done  by  modifying  the  r?2  parameter  as  a 
function  of  the  level  of  activity.  The  modeling  effort  focused  on  three  issues:  (I)  what  metric 
should  be  used  to  characterize  the  level  of  geomagnetic  activity  ( D,t  Kv  etc.),  (2)  what  is 
the  effect  of  changes  in  this  metric  in  the  occurrence  statistics,  and  (3)  how  can  this  effect 
be  modeled. 

In  the  work  to  determine  which  metric  to  use,  we  focused  on  the  following  parameters: 

1.  Kp  at  the  time  of  the  scintillation  observation, 

2.  Kp  at  sunset  (<e  =  0,  te  as  defined  in  SR4), 

3.  D,t  at  the  time  of  the  scintillation  observation, 

4.  Minimum  D,t  observed  during  the  previous  day  and  the  local  time  of  the  minimum. 


Although  D,t  is  not  currently  available  in  either  real-time  or  in  forecasts  at  the  USAF  Space 
Forecast  Center  (a  primary  customer  for  the  WBMOD  upgrade),  work  is  under  way  to 
obtain  this  parameter  for  use  in  other  space  models.  The  fourth  item  is  an  attempt  to  model 
the  behavior  described  in  Aarons  [1991]  in  which  the  local  time  that  a  storm  reaches  its 
minimum  D,t  controls  the  behavior  of  nighttime  scintillation.  After  some  work  with  these 
two  parameters  (the  minimum  D$t  and  time  of  the  minimum),  it  became  apparent  that  while 
they  could  possibly  be  useful  in  a  phenomenological  model,  the  changes  seen  in  the  occurrence 
statistics  when  categorized  using  these  parameters  were  not  that  much  different  (or  clearer) 
than  the  changes  seen  in  the  simpler  parameters  in  the  above  list.  When  a  successor  model 
to  WBMOD  that  deals  with  scintillation  from  a  phenomenological  standpoint  is  developed, 
this  set  of  parameters  should  be  revisited. 

Figures  25  through  27  are  samples  of  what  we  found  in  this  study.  These  figures 
show  the  occurrence  statistics  of  S4  >  0.5  as  a  function  of  te  and  day  of  the  year  from  the 
Huancayo  VHF  data  set  for  three  ranges  of  D,t  (Figure  25),  Kp  (Figure  26),  and  I\p  at 
sunset  (Figure  27).  (The  general  behavior  seen  in  the  Huancayo  data  is  fairly  representative 
of  all  three  stations.)  All  three  figures  show  the  same  general  behavior  as  the  level  of  activity 
increases: 
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Figure  27:  Occurrence  statistics  of  S4  >  0.5  from  the  Huancayo  VHF  data  set  as  a  function 
of  TPSS  (<e)  and  day  of  the  year  for  three  levels  of  geomagnetic  activity  as  measured  by 
A'pat  te  =  0:  (a)  Kp  <  2o,  (b)  2o  <  Kp  <  4o,  and  (c)  Kp  >  4°. 


1.  The  occurrence  of  scintillation  at  the  post-sunset  near-equinoctal  peaks  decreases  (re., 
inhibition  of  post-sunset  plume  development). 

2.  The  post  sunset  occurrence  statistics  in  the  local-summer  minimum  between  the  peaks 
do  not  change  as  much,  giving  rise  to  the  appearance  that  the  two  peaks  are  “merging” 
into  a  single  peak  in  local  summer. 

3.  The  variation  of  the  occurrence  statistics  with  time  become  “stretched  out”  in  the 
sense  that  the  levels  become  higher  in  the  post-midnight  period. 


In  reviewing  these  figures,  we  came  to  the  following  conclusions: 


1.  For  a  climatological  model  of  scintillation  occurrence  statistics,  the  l\p  parameter  is 
as  good  as  Dsi  (or,  at  least,  no  worse)  for  characterizing  the  effects  of  geomagnetic 
activity  on  eq-  atorial  scintillation. 

2.  The  post-sunset  reduction  in  the  occurrence  frequency  shows  most  clearly  in  the  I\p 
at  sunset  parameter,  and  the  post-midnight  enhancement  in  the  occurrence-frequency 
levels  shows  more  clearly  in  the  Kv  valid  at  the  time  of  the  observation. 


These  observations  will  form  the  basis  for  the  model  variation  of  the  log(C*/,)  PDF  with 
geomagnetic  activity. 

A  fmal  'ote  —  improvements  to  this  aspect  of  the  model  await  a  better  understanding 
of  the  exact  details  of  the  coupling  between  processes  at  work  in  the  magnetosphere  and 
at  high  latitudes  during  geomagnetic  storms  and  the  observed  behavior  of  scintillation  in 
the  equatorial  region.  Recent  papers  by  Aarons  [1991]  and  Kelly  and  Maniyainu  [1992]  are 
definite  steps  towards  a  better  understanding,  but  further  work  is  needed. 

4.  Propagation  Model  Development 

The  propagation  model  in  WBMOD  has  a  known  deficiency  in  the  way  in  which  S4  is 
calculated  for  a  two-way  path.  The  present  work  is  focused  on  calculating  the  expected 
correlation  between  scintillation  imposed  on  a  signal  on  the  up-  and  down-links  of  a  two-way 
path,  'flie  basic  equation  relating  S4  calculated  on  a  one-way  path  (denoted  Sm)  and  that 
for  the  two-way  path  is 


si  + 1 

pSl  +  1 


1  +  ps‘  4  +  2 


PSl 


SI  F  1 


(13) 


where  p  is  the  correlation  between  the,  up-  and  down-links  [ Fremouw  and  Ishimaru ,  1992]. 
In  the  initial  implementation  of  this  calculation,  we  plan  to  implement  Equation  [13]  with 
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a  user-supplied  value  for  p  ranging  from  0  (full  decorrelation)  to  1  (full  correlation).  1  he 
guidance  provided  to  the  user  will  be  to  enter  1  for  full  correlation  (which  will  be  true  in 
most  cases)  unless  there  is  good  reason  to  enter  some  value  less  than  1.  When  the  theory  for 
calculating  p  from  the  system  geometry  and  the  effective  scan  velocity  (K//)  is  completed, 
that  theory  will  be  implemented  in  an  upgrade  to  WBMOD  and  the  user-input  value  for  p 
will  serve  as  an  optional  user-specified  override  to  the  value  calculate*1  internally. 

5.  Conclusion 

The  next  stage  in  the  upgrade  of  the  equatorial  model  will  be  completion  of  modeling  the 
variation  of  scintillation  with  geomagnetic  activity  and  inclusion  of  the  variation  with  the 
solar  cycle  (i.e.,  with  SSN).  When  the  modeling  is  completed,  the  new  model  will  be  merged 
into  WBMOD,  and  the  user  interface  (inputs  and  outputs)  will  be  modified  as  necessary  for 
the  new  model  requirements  and  capabilities. 

Once  the  equatorial  modeling  is  completed,  the  effort  will  move  to  the  high  latitude 
region  using  data  from  the  Wideband,  HiLat,  and  Polar  BEAR  data  sets. 
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